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Abstract. We describe the main-sequence evolution of 
a rotating 9 M Q star. Its interior rotation profile is de- 
termined by the redistribution of angular momentum 
through the meridian circulation and through the shear 
turbulence generated by the differential rotation; the pos- 
sible effect of internal waves is neglected. We examine the 
mixing of chemicals produced by the same internal mo- 
tions. Our modelization is based on the set of equations es- 
tablished by Zahn (1992) and completed in Matias, Talon 
& Zahn (1996). Our calculations show that the amount 
of mixing associated with a typical rotation velocity of 
~ 100 kms -1 yields stellar models whose global parame- 
ters are very similar to those obtained with the moderate 
overshooting (d/Hp ~ 0.2) which has been invoked until 
now to fit the observations. Fast rotation (~ 300 kms -1 ) 
leads to significant changes of the C/N and O/N surface 
ratios, but the abundance of He is barely increased. The 
modifications of the internal composition profile due to 
such rotational mixing will certainly affect the post-main- 
sequence evolution. 

Key words: Stars: abundances, early-type, evolution, 
interior, rotation. 



1. Introduction 

There are many reasons to suspect that the amount of 
mixing presently included in stellar models, in the form 
of overshooting or semi-convection, is insufficient to ex- 
plain some observational facts (Maeder 1995). In partic- 
ular, Herrero et al. (1992) found a strong correlation be- 
tween He enrichment and surface velocity in O-type stars, 
which indicates that there is probably a link between ro- 
tation and mixing. 

Send offprint requests to: S. Talon 



Indeed it is well known, since the early work of Edding- 
ton (1925) and Vogt (1925), that the meridional currents 
arising in a rotating star may cause internal mixing. But 
until recently, the only case considered seriously was that 
of uniformly rotating stars, and it has been proved by Mes- 
tel (1953) that in such stars the stabilizing effect of the 
composition gradients, which arise from hydrogen burn- 
ing, is so strong that even the fastest rotators are unlikely 
to be mixed. 

The situation is different, however, if one allows the 
circulation to modify the rotation profile, as has been 
pointed out by Zahn (1992; Paper I). He invoked the 
anisotropic turbulence observed in stratified fluids, with 
much stronger transport in the horizontal than in the ver- 
tical direction, to reduce the problem to one dimension, 
with the angular velocity f2 depending only on depth. 
Then the evolution of VL obeys a "hyper-diffusion" equa- 
tion, involving the fourth derivative of f2, and in massive 
main-sequence stars the rotation profile tends to adjust 
such as to minimize the flux of angular momentum, as 
had been anticipated by Busse (1981). Provided there is 
no interference from other transport processes, the resid- 
ual turbulence associated with that differential rotation is 
the main cause of mixing in the radiative interior, and this 
is true even in the regions of strong composition gradient 
which arise near the convective core, thanks to the ho- 
mogenizing action of the anisotropic turbulence, as shown 
recently by Talon & Zahn (1996). 

Other groups have worked on rotational mixing in 
early-type stars, and they have clearly shown that such 
mixing may modify the surface abundances of He, C, N 
and O (Langer 1991, 1992; Denissenkov 1994; Eryurt et 
al. 1994). However, they have used an older prescription 
based largely on dimensional arguments, which links that 
mixing to the rotation velocity (Zahn 1983). Meynet & 
Maeder (1996) performed similar calculations, but they 
took into account differential rotation, and calculated the 
diffusion related to it, assumimg that the rotation profile 
is determined by the local conservation of angular mo- 



2 



S. Talon et al.: Rotational mixing in early-type stars 



mentum. The improvement brought by the present paper 
is that here the internal rotation evolves for the first time 
consistently under the combined action of the meridian 
circulation and the shear turbulence. Furthermore, as ex- 
plained in Zahn (1992) and in Matias et al. (1996), our 
treatment of rotational mixing involves only one free pa- 
rameter, which we hope to calibrate in the near future 
using fast rotators. 

In this paper our goal is more modest, for we shall not 
consider stars as massive as those observed by Herrero et 
al. The main reason is that we do not want our results 
to dependent too drastically on the prescription employed 
to treat semi-convection, which significantly affects the 
structure of the inner portion of O-stars. We choose to 
describe the evolution of a 9 M Q star, and intend to com- 
pare our results with the observations used to calibrate 
overshooting in moderate mass stars. 

We begin by presenting our stellar models and the 
equations governing the transport of angular momentum 
and of chemical elements (§2). The results for a hypothet- 
ical, stationary model are discussed in (§3), and those for 
the full main-sequence evolution in (§4). 

2. Stellar models with rotation 

2.1. Input physics 

Our models were built with the Geneva stellar evolution 
code, and they were evolved starting from a homogeneous 
composition at the zero-age main-sequence (ZAMS). 

The equation of state includes partial ionization close 
to the surface of the star. 

We used the OPAL radiative opacities (Iglesias et al. 
1992) which include the spin-orbit interactions for Fe. The 
relative metal abundances are based on Grevesse (1991), 
consistently with the opacity used. Note however that this 
consistency will be lost in the course of evolution, since 
mixing and nuclear processing modify those relative abun- 
dances. The tables are completed for the lower tempera- 
tures (below 10 000K) by those given by Kurucz (1991). 

Hydrogen burning proceeds through the p-p chains and 
the CNO cycle. The reaction rates are as given by Caugh- 
lan & Fowler (1988) with the correction by Landre et al. 
(1990) for the reactions 17 0(p,7) 18 F and 17 0(p,a) 14 N. 

For the mixing length parameter a — I /Hp we take 
the value 1.6, drawn from the solar calibration with the 
same code (cf. Schaller et al. 1992). 

No allowance will be made here for convective pene- 
tration, which is still a matter of debate. Most theoret- 
ical predictions (Roxburgh 1978, 1989; Zahn 1991) tend 
to overestimate this penetration, and the reason may be 
that they do not account for the inhibiting effect of rota- 
tion, which has been clearly demonstrated in the recent 
numerical simulations performed by Julicn et al. (1996). 

One important observational constraint for models of 
OB stars is the width of the main-sequence, and it can be 



accounted for by an overshooting parameter d/H p ~ 0.2 
(cf. Stothers & Chin 1991). However, this value charac- 
terizes in fact the size of the well mixed core of the star, 
and that mixing may be due either to convective penetra- 
tion or to any other mechanism. In this study, we shall 
examine whether this core extension may be attributed to 
rotational mixing, as described by Zahn's formalism (cf. 
Paper I). Therefore, we assume that convective overshoot- 
ing is rather weak, just enough to yield a finite value of 
the subadiabatic gradient (V a <a — V ra d), in order to prevent 
the divergence of the circulation velocity at the boundary 
of the convective core. 

Mass loss is also included, following the empirical re- 
lation given by de Jager et al. (1988), which is valid for O 
to M main-sequence stars. 

Furthermore, the equilibrium structure of our models 
includes the mean centrifugal force (averaged over an iso- 
bar). The hydrostatic equation reads: 
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where we use standard notations for the pressure (P), den- 
sity (p) , radius (r) , mass (m) , gravitational constant (G) , 
colatitude (9), angular velocity (f2) and where the hori- 
zontal average of / is defined as 
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We shall not refine here the implementation of the cen- 
trifugal force, since it would have little effect on mixing; for 
a more rigorous treatment, we refer the reader to Meynct 
& Maeder (1996), where they adapt the Kippenhahn & 
Thomas (1970) method to incorporate the hydrostatic ef- 
fects of rotation in one dimension to the case of shellular 
rotation. 

2.2. Transport equations 

In this paper, we shall consider only the transport of mat- 
ter and angular momentum due to meridional circulation 
and shear turbulence. Microscopic diffusion has negligible 
effect in early-type stars, and it will be ignored here. The 
possible redistribution of angular momentum through in- 
ternal waves is also neglected, although it has been shown 
to play an important role in the Sun (Kumar & Quatacrt 
1996; Zahn, Talon & Matias 1996); we shall revisit the 
problem once we dispose of a reliable prescription for this 
transport in a rapidly rotating star. 

The equations controlling the evolution of the rotation 
profile through meridional advection and turbulent diffu- 
sion have been presented in Paper I and in Matias et al. 
(1996). Let us recall that the main assumption in this de- 
scription is the strong anisotropy of the turbulence, whose 
horizontal transport enforces a state of shellular rotation 
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£1 = fi(.P), so that all variables may be expressed with 
respect to the isobars as 



f(P,6) = (f(P,6))+f(P)P 2 (co S 6) 

= /(p) + /(p)p 2 (cos e). 



(3) 



After performing the suitable horizontal means, the 
transport equation for angular momentum is cast into 
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(cf. Eq. (2.7) of Paper I), where v v is the vertical (turbu- 
lent) viscosity and U the amplitude of the vertical circu- 
lation speed u(P,6) = U{P)P 2 {cos6): 
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with L being the luminosity, g the local gravity, Cp the 
specific heat, T the temperature. The function En de- 
pends on the velocity profile 
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and Eft, on the chemical inhomogeneities along isobars 
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(cf. Matias et al. 1996). Here e is the nuclear energy gen- 
eration rate, Ht the temperature scale height, \ the ra- 
diative conductivity, and c, the concentration of a chemi- 
cal species of mass number i; the subscript m designates 
the mean value of the considered variable over the iso- 
bar, and other subscripts refer to logarithmic derivatives. 
6 = r 2 /igdVi 2 / dr is a measure of the differential rota- 
tion in the vertical direction and A = p/ p is the variation 
amplitude of the mean molecular weight along an isobar. 
This amplitude is derived from the horizontal variation of 
the concentration Cj of each species, which results from 
the competition between horizontal diffusion and vertical 
advection: 
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The consequence is that the vertical transport due to the 
meridian circulation is turned into a diffusion, with an 
effective diffusivity given by (Chaboyer & Zahn 1992): 

\rU{r)\ 2 



D eS = 
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30 D h 

Thus the transport of a given chemical element varies 
with time according to 
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where the nuclear production/destruction rate c, is taken 
from the stellar structure code. 

In this study, the only contribution retained in the ver- 
tical turbulent diffusivity (D v ) is that of the vertical shear, 
and we assume that this diffusivity equals the turbulent 
viscosity. We use the prescription given by Talon & Zahn 
(1996): 

8Ri c {rdn/drf 



D„ = 



5 N 2 /(K + D h ) + N 2 /D h ' 
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which takes into account the homogenizing effect of the 
horizontal diffusion (Dh) on the restoring force produced 
by the ^-gradient. Here, the Brunt- Vaisala, frequency has 
been split into iV| = (V ad - V) and A 2 = -g^V M , K 
is the radiative diffusivity and Ri c w 1/4 is the critical 
Richardson number. 

It remains to evaluate the magnitude of the horizontal 
diffusion (Dh), which we are unable to derive from first 
principles. We will use a parametric relation established 
as follows. Since the horizontal shear is sustained by the 
advection of momentum, Dh must be related to the circu- 
lation velocity U in order to keep the differential rotation 
in latitude below a certain level. This leads to 



D h 
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where Ch is an unknown parameter of order unity (see 
Zahn 1992). 

2.3. Numerical technique 

At each time step, the stellar evolution code is used to 
compute the mean structure of the model. Then, the 
meridional circulation is calculated, with the induced dif- 
fusion. The new rotation profile, as well as the new chem- 
ical abundances, are then implemented back into the evo- 
lutionary code, and we iterate until we reach the solution. 

The fourth order equation governing the evolution of 
the rotation profile is solved by a relaxation method based 
on finite differences, as in the Hcnycy method employed 
in the stellar evolution code (cf. Henyey et al. 1964). 

The diffusion equation is solved from the center of the 
star up to the surface, with a finite element method, using 
a semi-implicit temporal integration. The convective zone 
is homogenized by diffusion, using an arbitrary (large) co- 
efficient of D — 10 14 cm s _1 . 
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Fig. 1. Normalized rotation profiles in the asymptotic regime, 
for a 9Mq star. The solid line represents the star with a surface 
velocity of 20 kms - , the long dashed line the same with a 
velocity of 100 kms -1 , the short dashed line with a velocity 
of 300 kms" 1 and the dotted line, the fastest rotator with a 
surface velocity of 425 kms -1 . 

3. Asymptotic regime 

When the star does not lose angular momentum, which to 
first approximation is the case of massive stars, the inter- 
nal rotation rate tends to a profile which cancels the flux 
of angular momentum in the radiative interior. If the star 
was not evolving, the flow would settle into a stationary 
regime, in which the advection of angular momentum by 
the circulation exactly balances the diffusive transport: 

- < 13 » 

Although this state will never be achieved in a real star, 
it is worth considering it for the purpose of comparison 
with other results. 

3.1. Asymptotic circulation profile in realistic models 

We calculated the stationary regime obtained for a ZAMS 
star of 9-M©. For this, we solved the p.d.e. (fourth order in 
space, first order in time) governing the evolution of angu- 
lar momentum (Eq. [|), starting from solid body rotation, 
until we reached a stationary regime. Urpin et al. (1996) 
performed a similar study, for 20M Q models. However, 
they used a procedure slightly different from ours, search- 
ing directly the solutions of the third order o.d.e. (|l3|). We 
shall compare our results to theirs when possible. 

In Fig. 1, we show the rotation profiles obtained 
for four different surface velocities (20, 100, 300 and 



Fig. 2. Circulation velocity in the radial direction U{r), for 
the asymptotic regime. (Line styles are defined as in Fig. 1.) 

425 kms -1 ). Figure 2 displays the corresponding radial 
velocities of the meridian circulation. 

We observe that the normalized profile is independent 
of the rotation velocity in the inner portion of the star, 
whereas the situation is quite different close to the sur- 
face. The reason is that v v scales as fi 2 , and so does also 
Eft, through g/g and O, as long as the next order term 
— £l 2 /2irGp remains negligible compared to unity, which is 
the case in the deep interior (cf. Eq. 6). Then Eq. (|l^) is 
insensitive to the magnitude of Q, and it can only deter- 
mine the shape of the rotation profile. That higher order 
term, which we owe to Gratton (1945) and Opik (1951), 
dominates however near the surface, and it renders Eq 
nonlinear in f2 2 , thus modifying the rotation profile. For 
rapid rotation, the radial derivative of becomes positive 
near the surface, imposing a reversal of the diffusive flux 
and therefore also of the radial component of the circula- 
tion velocity, as can be seen in Fig. 2. 

Urpin et al. reach similar conclusions. They obtain a 
differential rotation between center and surface of ~ 1.15, 
which is comparable to what we get. Their fastest rota- 
tor (#4) corresponds to our model with 300 kms -1 . Its 
circulation velocity becomes marginally negative near the 
surface, as for ours. However, the velocities they obtain 
for the inner part of their models (<~ 50 cms -1 ) imply 
viscosities of order ~ 10 12 cm 2 s -1 , values that we do not 
reach. 

Figure 3 displays the diffusion coefficients D = D e ff + 
D v + -D ra d in the asymptotic regime (-D r ad is the radia- 
tive viscosity). In the 300 and the 425 kms -1 models, we 
notice a sharp drop in D around r/R ~ 0.95 and 0.85 
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Fig. 3. Diffusion coefficients in the asymptotic regime. (Line 
styles are defined as in Fig. 1.) 



Fig. 4. Rotation profile required by thermal equilibrium; com- 
parison between a ZAMS model and a TAMS model. 



respectively. It corresponds to the place where the deriva- 
tive of vanishes and thus, where no turbulence occurs. 
There, the viscosity is purely radiative and the transport 
of heavy elements will be strongly hindered. One could be 
tempted to conclude that, for the fastest stars, there will 
be two partially mixed zones, separated by a gap where 
very little diffusion takes place. Even worse, the region just 
outside of the convective core is not turbulent either, the 
diffusion coefficients diminishing by one or two orders of 
magnitude there. But we shall see in the next section that 
these conclusions are strongly affected by the structural 
changes accompanying the evolution, in particular by the 
chemical composition gradients. 

Figure 3 illustrates yet another property of such rotat- 
ing models, namely that the size of their convective core 
depends on rotation, with the fastest model having the 
smallest core. This is in agreement with previous studies 
of rotating stars (cf. e.g. Clement 1979), which predicted 
a lower central temperature and a smaller convective core 
for a rotating star than for a non-rotating one. 

4. Circulation and evolution 

The problem described in the previous section is rather 
academic. It is merely a guide to understand some basic 
properties of the circulation. However, it gives no clue as 
to what will happen actually to an evolving model, and it 
bypasses the important role of the /x-gradients. 

As is well known, the radius of early-type stars in- 
creases as they evolve on the main-sequence. However, 
this expansion does not proceed in an homologous way, 



for the inner part is shrinking instead. Local conservation 
of angular momentum would lead to a core rotating much 
faster than the surface. Since the profile of £1 required by 
the asymptotic regime does not change much between con- 
secutive models, we thus expect a circulation to arise that 
transports momentum from the center of the star to the 
surface. Such a circulation will rise along the equator and 
descend along the pole. Note that in a rigidly rotating star 
the circulation would be in the opposite direction (Sweet 
1950). 

In Fig. 4 we show the rotation profile required for this 
hypothetical stationary regime at the end of the main- 
sequence, compared to that on the ZAMS, for a model 
with an initial surface velocity of 100 kms -1 . It illustrates 
the need to carry angular momentum outward, since local 
conservation of angular momentum would lead to a dif- 
ferential rotation ratio (f2 c /f2 s ) of order ~ 10. (Note that 
on the TAMS, the total radius of the star has doubled 
compared to the ZAMS model). 

Mass loss is also an important cause of circulation. The 
constant readjustment of the surface layers, which have to 
keep moving outward as the star loses mass, involves an 
important circulation in the outer part of the envelope. 
We thus expect outer layers to remain well mixed. 

As was mentioned before though, the main impact of 
evolution is felt through the formation of a mean molec- 
ular weight gradient. In an inhomogeneous star, the rele- 
vant quantity in the description of the circulation is not 
ft itself, but rather 9 — A (cf. Eq. 5). As was first shown 
by Mestel (1953), the advection of matter of larger mean 
molecular weight from the center of the star has a chok- 
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Fig. 5. The evolution of a rotating star of 9 Mq in the HR 
diagram. The solid line tracks the evolution of a star with an 
initial surface velocity of 100 kms -1 , the long dashed line that 
of a star with an initial surface velocity of 300 kms -1 . For 
reference, the short dashes show the evolution of the same star 
with an overshooting of d/H p = 0.2 and the dotted line the 
star without any overshooting. 

ing effect on the circulation. This, of course, is particu- 
larly important near the regressing convective core. How- 
ever, contrary to Mestel's conclusions, the circulation is 
not completely "killed" by those gradients, even though 
we expect it to be slowed. In an evolving star, the hori- 
zontal //-gradients keep adjusting themselves to the new 
equilibrium state; furthermore, due to the erosion by the 
horizontal turbulence, at least a little amount of advection 
is required to maintain the horizontal inhomogeneities. 

Another critical effect of the /x-gradients is that they 
tend to inhibit baroclinic shear instabilities. However, as 
was shown by Talon & Zahn (1996), if the turbulence 
is strongly anisotropic, horizontal diffusion reduces the 
buoyancy force associated with these gradients, just like 
temperature gradients are weakened by thermal diffusion. 
Their prescription ( |TT| ) allows for vertical mixing, provided 
the differential rotation is strong enough. 

4--1- Evolutionary calculations 

Complete evolutionary calculations were performed for 
two models, initially rotating at 100 and 300 kms -1 (we 
will later refer to these models as "slow" and "fast"). 
The stars are assumed to be in the asymptotic regime 
on the ZAMS, thus having the rotation profiles shown in 
Fig. 1. We are aware that this initial condition is arbitrary, 



Fig. 6. Evolution of the central abundance of hydrogen. (Line 
styles are defined as in Fig. 5.) 

but the memory of the initial condition will be lost in an 
Eddington-Sweet time-scale 

GM 2 ( GM \ . A . 

tES = ^{wW)' (14) 

of order 10 6 years for the considered models. Furthermore, 
during this early adjustment phase, no chemical inhomo- 
geneities are present, and this diminishes the impact of 
the initial condition on the amount of mixing. 

In Fig. 5, we show the evolutionary paths followed by 
the two models, plus that of two comparison models, one 
with a convective overshooting of d/H p = 0.2, and the 
other one without overshooting. One striking feature of 
that diagram is that the mixing associated with a star ini- 
tially rotating at 100 km s -1 yields an evolutionary path 
which can hardly be distinguished from that of a model 
with a small amount of overshooting. The initial positions 
of these models on the HR diagram also illustrate the small 
impact of the hydrostatic corrections on the global param- 
eters, compared to those related to mixing. 

At the end of the main-sequence, the "slow" model has 
lost 0.09% of its initial mass, and 1.4% of its total angular 
momentum. The "fast" model has lost 0.12% and 2.3% of 
its mass and momentum respectively. Remember that we 
used an empirical formula for mass loss that does not take 
into account the rotation rate of the star. The difference in 
the mass lost only reflects different lifetimes and effective 
temperatures during the evolution. 

Figure 6 displays the evolution of the central mass frac- 
tion of hydrogen, and thus also the lifetime on the main- 
sequence. For the "slow" model, which has a velocity typ- 
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Fig. 7. Total diffusion coefficient and its variation with evolu- 
tion, (a) Model with an initial surface velocity of 100 km s~ . 

(b) Model with an initial surface velocity of 300 km s" 1 . The Fi S- 8 - Rotation profile and its variation with evolution. (Line 
full line refers to the youngest model, the dotted line to the st y les are defined as in Fig. 7 and refer to the same models), 
model in the middle of the main- sequence, and the dashed line 
to the oldest model, at the turning point of the main-sequence. 



ical for observed 9 Mq stars, the evolution in the HR dia- 
gram, as well as the main-sequence lifetime, are compara- 
ble to what is obtained with a standard model including 
overshooting. 

As mentioned in the Introduction, convective penetra- 
tion may be severely hindered in rotating stars, accord- 
ing to the recent simulations carried out by Julien et al. 
(1996). Our calculations show that the widening of the 
main-sequence may be due to rotational mixing, in stars 
having a typical velocity (~ 100 km s _1 ). 

4-2. Role of fi- gradients 

In the preceding section, we described the global results 
obtained by solving the complete set of equations stated in 
§2, which include the effect of the composition gradients. 
Here we shall examine in more detail the particular role 
of these /i-gradicnts. 

In Fig. 7, we show the diffusion coefficient at three dif- 
ferent ages, for the two rotating models. For the "slow" 
model, the major contribution to diffusion is turbulent. 
Close to the convective core, we note a significant drop in 
the diffusion coefficient, which is due to the steep molec- 
ular weight gradient. The turbulent transport is affected, 
but it is not suppressed (cf . Eq. [ll]) . 



< 
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Fig. 9. A = and its variation with evolution. (Line styles 
are defined as in Fig. 7 and refer to the same models). 
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Fig. 10. Evolution of surface abundances. (Line styles are de- 
fined as in Fig. 5). 

For the "fast" model, we note also a strong diminu- 
tion of the diffusion coefficient near the surface. It is due 
to the reversed circulation cell, which grows with time. 
The reason for this growth is readily found; as the star 
evolves off the zero-age main-sequence, its radius grows 
and its breakup velocity decreases. Local conservation of 
angular momentum would lead to still faster decline of 
the surface velocity. However, circulation and turbulence 
keep extracting angular momentum from the interior, and 
the actual ratio between surface and critical velocities de- 
creases. 

The minimum of the total diffusion coefficient is thus 
located where dil/dr = 0, and it reduces there to the 
effective diffusivity D e g associated with the meridian cir- 



Fig. 11. Final abundance profiles on the TAMS of hydrogen 
(H), helium (He 4 ), carbon (C 12 ), nitrogen (N 14 ) and oxygen 
(O 16 ). (Line styles are defined as in Fig. 5). 



culation (cf. ||) . That strangling will be the limiting factor 
in the modification of the surface abundances of the fast 
models (see also Figs. 10 and 11). 

Another important effect is shown on the rotation pro- 
file within the model (cf. Fig. 8). It is clear that those pro- 
files are quite different from what is required for thermal 
equilibrium (see Fig. 4). The main reason for this differ- 
ence is the presence of the term in A (shown in Fig. 9). 
The steep profile in A requires that be also very steep 
there, so that the sum of the two terms can give a profile 
similar to that shown in Fig. 4, which would be required 
for the asymptotic regime. 
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There is here an interesting coupling between all the 
variables. If d/z/dr is large, then so is A (cf. Eq. and 
consequently, so is dlnfi/dr. Such a strong differential 
rotation enhances turbulent diffusion, thereby diminishing 
d/i/dr. However, as turbulent diffusion increases, so does 
U r , and thus so does also A. 

Let us stress that the behavior of A depends on an 
adjustable parameter Ch of order unity. That parameter 
is also present in D v , through Dh (cf. Eq. [ll]). It is in- 
teresting to note that as Ch increases, so does A, but D t 
varies in the opposite way. The net result is that the final 
profiles depend sensitively on the choice of Ch, which has 
here the value 0.15. 

In Fig. 10, we show the evolution of the surface abun- 
dances of a few elements, and in Fig. 11, the final abun- 
dance profiles. Fast rotation is required to yield detec- 
tectable abundance changes on the main-sequence. How- 
ever, even for the "slow" model, the post-main-sequence 
evolution is probably affected by the composition changes 
inside the star. Furthermore, as mass loss becomes more 
important with the evolution off the main-sequence, ad- 
ditional changes of the surface abundances will occur 
through the "peeling" of the outer layers. 

5. Discussion 

We completed a first set of calculations of rotational mix- 
ing induced by meridional circulation as well as shear 
turbulence, with the rotation profile evolving consistently 
with these flows. 

Unfortunately, very few observations permit to verify 
the full validity of the results. Gies & Lambert (1992) have 
observed galactic B stars and measured their rotational 
velocity, as well as their CNO abundances. However, they 
studied only slow rotators (v sini < 85kms _1 ), which 
did not show strong evidence of composition changes. 

One would expect binary systems to yield important 
information about the precise position in the HR dia- 
gram of stars with known masses, radii and rotation ve- 
locities. The best known binary in the 9 M Q mass range 
is QX Car, whose masses are 9.267 ± 0.122 Mq and 
8.480 ± 0.122 Mq (cf. Andersen 1991). However, these 
stars are located too close to the ZAMS to give any inter- 
esting information about rotational mixing. Other binary 
systems (e.g. a Vir, NY Cep) are not known with sufficient 
precision to be used for any reliable comparison. 

The most promising tests will be provided by the 
fastest rotating B-stars, which show signs of helium en- 
richment. The difficulty will be however to sort out the 
effects of semi-convection from those of rotational mixing. 

Let us finally recall that a complete treatment of rota- 
tional mixing must include the angular momentum trans- 
port by the internal gravity waves emitted from the con- 
vective core. Such transport appears to be quite important 
in the Sun, as pointed out recently by Kumar & Quataert 
(1996) and by Zahn, Talon & Matias (1996). Since inter- 



nal waves tend to restore uniform rotation, which would 
speed up the circulation, one expects a priori more mix- 
ing to occur. But the whole problem is highly non-linear, 
and it is difficult to make any prediction until this wave 
transport has been implemented and its impact has been 
evaluated. Work is progressing in this direction. 
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